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Abstract 

Many of the factors that contribute to CD8 + T cell immunodominance hierarchies during viral infection are known. However, 
the functional differences that exist between dominant and subdominant epitope-specific CD8 + T cells remain poorly 
understood. In this study, we characterized the phenotypic and functional differences between dominant and subdominant 
simian immunodeficiency virus (SIV) epitope-specific CD8 + T cells restricted by the major histocompatibility complex (MHC) 
class I allele Mamu-A*01 during acute and chronic SIV infection. Whole genome expression analyses during acute infection 
revealed that dominant SIV epitope-specific CD8 + T cells had a gene expression profile consistent with greater maturity and 
higher cytotoxic potential than subdominant epitope-specific CD8 + T cells. Flow-cytometric measurements of protein 
expression and anti-viral functionality during chronic infection confirmed these phenotypic and functional differences. 
Expression analyses of exhaustion-associated genes indicated that LAG-3 and CTLA-4 were more highly expressed in the 
dominant epitope-specific cells during acute SIV infection. Interestingly, only LAG-3 expression remained high during 
chronic infection in dominant epitope-specific cells. We also explored the binding interaction between peptide:MHC 
(pMHC) complexes and their cognate TCRs to determine their role in the establishment of immunodominance hierarchies. 
We found that epitope dominance was associated with higher TCR:pMHC affinity. These studies demonstrate that 
significant functional differences exist between dominant and subdominant epitope-specific CD8 + T cells within MHC- 
restricted immunodominance hierarchies and suggest that TCR:pMHC affinity may play an important role in determining 
the frequency and functionality of these cell populations. These findings advance our understanding of the regulation of T 
cell immunodominance and will aid HIV vaccine design. 
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Introduction 

Virus-specific CD8 + T cells contribute to the control of HIV 
and SIV replication [1,2] and are therefore an important 
component of a protective immunity. In infected or vaccinated 
individuals, the frequencies of different viral epitope-specific CD8 + 
T cells vary considerably [3,4]. Epitopes restricted by the same 
MHC Class I allele can be ranked in an immunodominance 
hierarchy based on the relative frequencies of their respective 
epitope-specific CD8 + T cells [5] . The determinants of immuno- 
dominance hierarchies have been explored in the past in efforts to 
enhance the magnitude of particular epitope-specific CD8 + T cell 
responses through vaccination [6,7]. However, such investigations 
have primarily focused on determining the mechanisms underly- 
ing the establishment of immunodominance [8- 12]; less is known 



about the resulting functional differences between dominant and 
subdominant epitope-specific CD8 + T cells. 

CD8 + T cell polyfunctionality is often associated with superior 
viral control [13]. Numerous studies have reported that individuals 
with superior control of HIV replication and delayed disease 
progression have higher frequencies of polyfunctional CD8 T 
cells and therefore such cells are thought to be an important 
component of a protective immune response during HIV infection 
[14-17]. Therefore, HIV vaccine development is particularly 
focused on the identification of vaccine strategies that will generate 
highly polyfunctional, and therefore protective, responses [18]. In 
order to do so, a better understanding of the factors that contribute 
to the polyfunctionality of a CD8 + T cell response is needed. 

Few studies have evaluated the functional differences between 
dominant and subdominant epitope-specific CD8 + T cells. Studies 
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Author Summary 

MHC-restricted CD8 + T cell populations that bind viral 
proteins are often present at different frequencies. It is 
thought that those virus-specific CD8 + T cells that are 
present at the highest frequency are predominantly 
responsible for eliciting control of viral infections. While 
the number of virus-specific CD8 + T cells is undoubtedly 
important, the functionality of these cells may also play an 
anti-viral role. It is not known if high-frequency virus- 
specific CD8 + T cells are more functionally effective against 
viral infection than those present at low frequencies. In this 
study, we characterized the functional differences between 
the SIV-specific cells present at high versus low frequencies 
in rhesus monkeys infected with simian immunodeficiency 
virus (SIV). We found that the high- and low-frequency SIV- 
specific cells had different functional capacities during 
acute and chronic SIV infection. We also found that the 
affinity with which a cell interacts with viral proteins may 
contribute to these functional differences. These findings 
further our understanding of anti-viral immune responses 
and may help to inform HIV vaccine development. 

in DNA-immunized mice showed that subdominant epitope- 
specific CD8 + T cells were less cytotoxic, but produced more 
cytokines, than dominant epitope-specific CD8 T cells during 
acute infection [19]. Moreover, a cryptic epitope-specific CD8 + T 
cell population exhibited an altered maturation phenotype when 
compared to cells comprising the dominant response [20] . A study 
of the functionality of GD8 + T cells within the HLA-B27-restricted 
immunodominance hierarchy in HIV-infected individuals found 
higher-frequency epitope-specific CD8 T cells were associated 
with superior in vitro viral suppression; although, subsequent 
functional studies could not identify a particular function 
associated with this suppressive capacity [21]. Additional knowl- 
edge of functional differences between dominant and subdominant 
epitope-specific CD8 + T cell responses is needed for enhancing 
CD8 + T cell epitope-specific responses by vaccination. 

In the present studies, we performed longitudinal gene 
expression analysis in SIV- infected Mamu-A*01 rhesus monkeys 
to assess potential functional differences between dominant and 
subdominant epitope-specific CD8 + T cells during acute infection. 
We demonstrate that dominant epitope-specific CD8 + T cells 
display gene expression patterns consistent with a more mature 
phenotype and harbor greater cytotoxic potential than subdom- 
inant epitope-specific CD8 + T cells during acute infection. Protein 
expression analyses of CDS T cells sampled during chronic 
infection confirmed the presence of these functional differences 
between dominant and subdominant epitope-specific CD8 + T 
cells. Finally, dominant and subdominant epitope TCR:pMHC 
affinities were correlated with relative immunodominance. These 
findings advance our understanding of the basis underlying 
immunodominance and may better inform T cell-based HIV 
vaccine design. 

Results 

pi 1 C- and p54AS-specific CD8 + T cell frequencies and 
viral loads during acute SIVmac251 infection 

We infected six Mamu-A*01 + rhesus monkeys with SIVmac251. 
The Mamu-A*0 1 -restricted dominant SIV Gag pllC- and 
subdominant SIV Env p54AS-specific CD8 + T cells were 
quantified and sorted by flow cytometry on days 14, 21, 28, 35, 
42, 56, and 70 post-inoculation. The Mamu-A*0 1 -restricted CD8 + 



T cell response is one of the best characterized CD8 + T cell 
responses in SIV-infected rhesus monkeys [22,23]. The pllC 
epitope (Gag CM9) has been used widely as a model dominant 
epitope in studies of SIV-specific immunity [24-26] . p54AS (Env 
TL9) is a well-described subdominant epitope that, from our 
experience, was the most likely to be elicited in all Mamu-A*01 
rhesus monkeys at a frequency that was high enough such that we 
could collect sufficient numbers of cells to facilitate our gene 
expression measurements. 

The differences in frequency between the pllC- and p54AS- 
specific CD8 + T cells were apparent upon detection in the 
peripheral blood (day 14) and were maintained throughout the 
duration of infection (Figure 1A and Figure SI). The plasma viral 
loads peaked at approximately 7 logs of viral RNA copies/mL by 
day 14 and reached setpoint of approximately 5.5 logs of viral 
RNA/mL by day 28 (Figure IB). 

We measured RNA expression on the sorted epitope-specific 
cell populations at various timepoints post infection. We also 
sorted total naive CD8 + T cells (CD95"CD28 + ) before chaUenge 
(day 0) and measured RNA expression in these cells to establish a 
baseline expression level for each transcript. We generated 
expression data from both epitope-specific cell populations from 
all six monkeys only on days 14, 21, 56 and 70. We obtained 
expression data for both epitope-specific populations in only some 
of the monkeys on the three other time points: days 35 and 42 (4 of 
6 monkeys) and day 28 (5 of 6 monkeys). Therefore, sometimes we 
observed trends suggesting differential expression (Figure S2) that 
did not achieve statistical significance possibly due to reduced 
power. 

To characterize the functional differences between the domi- 
nant pllC- and subdominant p54AS-specific cells based on their 
RNA expression profiles, we focused expression analysis on genes 
known to be associated with specific CD8 + T cell functions. These 
genes were grouped into the following categories: maturation, 
cytotoxicity, cell cycle and apoptosis, and cytokines and chemo- 
kines (Table 1). Genes were considered differentially expressed if 
they satisfied both of the following criteria: 1) expression was 
significandy different between the two epitope-specific cells on at 
least one timepoint and 2) the median fold difference of expression 
at this timepoint (pllC/p54AS) was at least ±1.5. Technical 
replicates were not conducted; gene expression was measured once 
per monkey per timepoint. 

Maturation-associated gene expression 

Of the genes involved in CD8 + T cell maturation, IL7R 
(CD127/IL-7 receptor, refseq accession # XM_937367.1), SELL 
(CD62L/L-selectin, NM_000655.3), and CCR7 (NM_00 1838.2) 
met our criteria for differential expression following SIV infection 
(Figure 2A). While the trends of expression of these three genes 
between pllC- and the p54AS-specific cells were similar over 
time, all of these genes were found to be expressed at relatively 
higher levels in the subdominant p54AS-specific cells compared 
to the dominant pi lC-specific cells. LL7R expression was higher 
in the p54AS-specific cells with significant fold differences of 1.8 
on day 14, 1.5 on day 56, and 1.9 on day 70. SELL expression 
was also higher in the p54AS-specific cells with significant fold 
differences of 1.8 on day 14, 2.0 on day 21, and 1.7 on day 56. 
CCR7 was more highly expressed in the p54AS-specific cells 
with a significant fold difference of 2.4 on day 14. While the 
expression of LL7R and SELL tended to be higher in the p54AS- 
specific cells for all timepoints examined, the expression for CCR7 
appeared only to be different on day 14 (Figure 2A and Figure 
S2A). 



PLOS Pathogens | www.plospathogens.org 



2 



April 2014 | Volume 10 | Issue 4 | e1 004069 



TCR Affinity and CTL Immunodominance Hierarchies 



A Tetramer Responses B Plasma Viral Loads 




Days Post Inoculation Days Post Inoculation 



Figure 1. Frequencies of pllC- and p54AS-specific CD8 T cells and plasma viral loads during primary infection. A) Mean frequencies 
of the p1 1 C- and p54AS-specific CD8 + T cells in peripheral blood. * significant at p<0.05. B) Plasma SIV RNA levels in the peripheral blood. 
doi:1 0.1 371 /journal, ppat.1 004069.g001 



Cytotoxicity-associated gene expression 

Of the genes involved in cytotoxicity, G^AIB (granzyme B, 
NM.004131.3), GZMK (granzyme K, NM_002104.2), and 
LAMP2 (CD107b, NM_013995.1) met our criteria for differential 
expression following SIV infection (Figure 2B). LAA1P2 was 
measured on the BeadChip by three different probes, two of 
which measured differential expression (probes 1659753 and 
1752351). Expression levels of each of the three cytotoxicity- 
associated genes were found to be increased in the dominant 
pi lC-specific cells compared to the subdominant p54AS-specific 
cells. On day 14, the pi lC-specific cells exhibited higher 
expression of both G%MB (significant fold difference of 1 .5) and 
G^MK (significant fold difference of 1.9). Although the trends in 
differential expression of these two granzyme genes were still 
apparent on day 21 and 28, differences were subsequentiy lost 
(Fig. 2B and Figure S2B). The two LAMP2 probes detected similar 
trends of expression of this gene. However, probe 1659753 
detected significant differential expression on day 14 (1.5-fold) 
while probe 1752351 detected significant differential expression on 
day 70 (1.8-fold). 

Cell cycle and apoptosis-associated gene expression 

Of the cell cycle and apoptosis-associated genes, nine met our 
criteria for differential expression following SIV infection 
(Figure 2C and Figure S2C). Two of these genes, BIRC3 (AIP1/ 
API2/c-IAP2, NM_001 165.3) and CASP2 (caspase-2, 
NM_032983.2), are involved in apoptosis. BIRC3, which encodes 
a protein with multiple anti-apoptotic functions [27], was 
differentially expressed on day 14, showing 1.8-fold higher 
expression in the p54AS-specific cells. CASP2, which encodes a 
pro-apoptotic protease [28], was differentially expressed on day 
14, showing 1.7-fold higher expression in the pi lC-specific cells. 

The remaining seven differentially expressed genes are involved 
in proliferation: UBE2C (ubiquitin-conjugating enzyme E2C, 
NM_181800.1),J£/it£B(auroraB kinase, NM_004217.2), CCND2 
(cyclin D2, NM_00 1759.2), CDC42 (cell division cycle 42, 
NM_001039802.1), MUSAP1 (nucleolar and spindle associated 
protein 1, NM_0 18454.5), UHRF1 (ubiquitin-like with PHD and 
ring finger domains 1, NM_00 1048201. 1), and ATM (ataxia 
telangiectasia, NM_000051.3). UBE2C and AURKB encode pro- 
proliferative genes involved in the regulation of the anaphase- 
promoting complex [29,30]; both were more highly expressed in 



the pi lC-specific cells on day 56 with 1.6- and 1.5- fold 
differences, respectively. CCND2, CDC42, MJSAP1, and UHRF1 
also encode pro-proliferative genes [31-35]. All of these genes 
showed differential expression on day 70 and were more highly 
expressed in the pi lC-specific cells: CCND2 by 1.6-fold, CDC42 by 
1.9-fold, mSAPl by 2.7-fold, and UHRF1 by 2.5-fold. Finally, 
ATM, an anti-proliferative serine/threonine protein kinase [36], 
was differentially expressed on day 21 and was more highly 
expressed in the p54AS-specific cells with by 1.5-fold. 

In summary, we detected differential expression of genes 
involved in CD8 + T cell maturation, cytotoxicity, and cell cycle 
and apoptosis suggesting that functional differences may exist 
between the dominant pi lC-specific and the subdominant 
p54AS-specific cells. Decreased expression of the genes encoding 
CD127, CD62L, and CCR7 and the increased expression of genes 
encoding granzyme B, granzyme K, and CD 107b (LAMP2) in the 
dominant pi lC-specific cells compared to the subdominant 
p54AS-specific cells suggested that the pi lC-specific cells were 
more mature and possessed greater cytotoxic potential than the 
p54AS-specific cells. Moreover, we observed that the pllC- 
specific cells had higher expression of pro-proliferative genes and 
an increased expression of a pro-apoptotic molecule. 

The dominant epitope-specific CD8 + T cell population 
was enriched with cells that exhibit a more mature 
phenotype 

Functional assays were performed on PBMCs isolated from 
chronically infected monkeys to confirm and extend the observa- 
tions of differential gene expression patterns. In addition to the 
dominant pllC- and subdominant p54AS-specific cells, we 
studied an additional CD8 + T cell population of even lower 
frequency that recognizes the Mamu-A*0 1 -restricted SIV Pol 
p68A (STPPLVRLV, Pol SV9 621 ) subdominant epitope [23]. 
Additionally, we analyzed cells from monkeys infected with a 
different SIV isolate, SIVsmE660 [37]. The p54 epitope from the 
SIVsmE660 (TVPWPNETL) virus differs by two amino acids 
from the p54 epitope from SIVmac251 (TVPWPNASL) virus [22] 
and we refer to this epitope as p54E660. The immunodominance 
hierarchy of the pllC-, p54AS/E660-, and p68A-specific cells 
during chronic SIVsmE660 infection was similar to SIVmac251- 
infected monkeys (Figure 3A and Figure S3). 
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Table 1. Genes compared between dominant and subdominant epitope-specific CD8 + T cells. 
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Figure 2. Genes differentially expressed between dominant pi 1C- and subdominant p54AS-specific CD8 + T cells. The median of the 
normalized expression values, measured in fluorescence units, of A) maturation, B) cytotoxicity, and C) proliferation and apoptosis genes that were 
found to be differentially expressed between the p1 1C- and p54AS-specific cells on at least one time point are shown for each timepoint evaluated 
during the first 70 days following SIVmac251 infection. The expression of these genes by total naive CD8 + T cells, measured on day 0, is also shown. 
For those genes with more than one probe on the BeadChip, the numerical probe IDs are included in the gene name. * indicates that expression met 
the criteria for differential expression on that timepoint. 
doi:1 0.1 371 /journal.ppat.1 004069.g002 



The maturation of tetramer-positive cell populations from 
SIVmac251- and SIVsmE660-infected monkeys was determined 
by measuring the surface expression of CCR7 and CD28. Cells 
were categorized as central memory (CCR7 + CD28 + ), transitional 
memory (CCR7 CD28 ), or fully differentiated effector memory 
(CCR7 CD28 ) as previously described for rhesus monkeys [38]. 
All tetramer-positive cells were CD95 + (not shown). The dominant 
pi lC-specific cells displayed a greater proportion of the more 
mature effector memory phenotype (CCR7 CD28 ) compared to 
the subdominant p54AS-.specific cells (Figure 3B and Figure S4). 
This is consistent with gene expression studies in SIVmac251- 
infected monkeys that expressed CCR7, CD62L, and SELL at lower 
levels in the pi lC-specific cells, suggesting they were more mature 
than the p54AS-specific cells. 

Additionally, we found that the more subdominant p68A- 
specific cells were phenotypically less mature compared to the 
subdominant p54AS/E660-specific cells (Figure 3B). Of the three 
epitope-specific cell populations, the p68A-specific population was 
the most enriched with cells of the less mature phenotypes and 
were the least enriched with cells of the mature effector memory 
phenotype. This trend was observed in both SIVmac251- and 
SIVsmE660-infected monkeys. While the distribution of pheno- 
types was similar between cells from SIVmac251- and 
SIVsmE660-infected monkeys, the sample size for SIVmac251 
(n = 4) was much smaller than that of SIVsmE660 (n = 8) and was 
not sufficiently powered to detect significance. 

To confirm that differences in gene expression of cytolytic 
molecules were similarly reflected in protein production, we 
measured the intracellular content of perforin and granzyme B 
within each tetramer-positive population. The dominant pllC- 
specific populations had a greater frequency of cells containing 
perforin and granzyme B than the subdominant p54AS/E660- 
specific population, which in turn had a greater frequency of 
perforin- and granzyme B-containing cells than the subdominant 
p68A-specific population (Figure 3C and Figure S5). The pllC- 
specific cells also contained more perforin and granzyme B per-cell 
(reflected by higher geometric mean fluorescence (GMF)) than the 
p54AS/E660-specific cells, which in turn contained more of these 
proteins than the p68A-specific cells (Figure 3D and Figure S5). 
These observations were consistent with previous findings that 
demonstrated not only that more mature cells are more likely to 
express perforin and granzymes [39], but also that more mature 
cells have more perforin and granzyme per cell than less mature 
cells [40]. 

The dominant p1 1C-specific CD8 + T cells exhibited 
decreased antigen-specific expansion compared to 
subdominant epitope-specific CD8 + T cells 

Although differences in expression of genes involved in cell cycle 
and apoptosis were observed, it was unclear what impact these 
genes had on cell expansion upon antigen stimulation during 
chronic infection. To explore this, PBMCs from monkeys chron- 
ically infected with SIVmac251 or SIVsmE660 were peptide- 
stimulated and expansion of the epitope-specific cells was assessed 
over 14 days. We found that the dominant pi lC-specific cells had 
reduced capacity to expand when compared to both of the 



subdominant epitope-specific cells (Figure 4 and Figure S6). This 
trend was consistent between cells from SIVsmE660- and 
SIVmac251 -infected monkeys and is in agreement with prior 
reports that effector memory CD8 + T cells have limited proliferative 
capacity compared to less mature memory CD8 + T cells [41]. 

The dominant p1 1C-specific CD8 + T cell population 
contained a lower frequency of cytokine- and 
chemokine-producing cells than the subdominant 
epitope-specific CD8 + T cell populations 

Our gene expression data did not demonstrate differences in 
cytokine or chemokine expression between the pi 1C- and p54AS- 
specific cells. However, cytokine genes often require re-stimulation 
immediately prior to measurement [42] . Therefore, we considered 
that differences in cytokine expression between dominant and 
subdominant epitope-specific cells might be observed following 
antigen stimulation. To investigate, we stimulated PBMCs with 
epitope peptides and measured IL-2, TNFa, IFNy, and MIPl-p 1 
production in standard intracellular cytokine staining (ICS) assays 
(Figure S7). Lower frequencies of cells expressing each of these 
cytokines and chemokines were present in the dominant pllC- 
specific cell population compared to the subdominant epitope- 
specific cell populations (Figure 5A). Furthermore, the subdomi- 
nant epitope-specific populations were more polyfunctional, 
exhibiting larger proportions of cells expressing three or four 
functions, compared to the dominant pi lC-specific population 
(Figure 5B and 5C). 

Expression of exhaustion-associated genes in dominant 
and subdominant CD8 + T cells 

To assess differences in functionality that may be due to 
different extents of exhaustion, we compared the expression of 
genes encoding exhaustion-associated molecules between the 
pllC- and p54AS-specific cells during acute infection (day 14), 
at the onset of chronic infection (day 70), and during late chronic 
infection (week 31 or 32). The genes LAG3 (NM_002286.4) and 
CTLA4 (NM_005214.3) were significantly more highly expressed 
in the dominant pi lC-specific cells during acute infection by 3.2- 
and 1 .5-fold, respectively (Table 2). However, no genes were found 
differentially expressed on day 70. During late chronic infection, 
LAG3 was expressed at increased levels in the dominant pllC- 
specific cells (1.7-fold). Therefore, the dominant pi lC-specific cells 
showed higher expression of exhaustion-associated genes than the 
subdominant p54AS-specific cells. 

High TCR:pMHC affinity correlates with the epitope 
dominance 

We were interested in identifying mechanisms underlying the 
differences in the frequency and function of epitope-specific cells 
within an immunodominance hierarchy. Since TCR:pMHC 
interactions can influence T cell activation, proliferation, and 
function [43-49], we considered that TCR:pMHC interactions 
might play a role in the establishment of the magnitude or of the 
functionality of CD8 + T cells within an immunodominance 
hierarchy. We employed a surface plasmon resonance (SPR)-based 
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Figure 3. Phenotype and cytotoxic potential of dominant and subdominant epitope-specific cells during chronic SIVmac251 and 
SIVsmE660 infection. The frequency, cell surface phenotype, and ex vivo perforin and granzyme B content were evaluated for pi 1C-, p54AS/E660-, 
and p68A-specific cells, identified with Mamu-A*01 tetramers, from peripheral blood of chronically infected SIVmac251- and SIVsmE660- infected 
rhesus monkeys. A) Frequencies of the pi 1C-, p54AS/E660-, and p68A-specific cells from (top) SIVmac251- and (bottom) SIVsmE660- infected rhesus 



PLOS Pathogens | www.plospathogens.org 



7 



April 2014 | Volume 10 | Issue 4 | e1 004069 



TCR Affinity and CTL Immunodominance Hierarchies 



monkeys. Error bars indicate median ± interquartile range. SIVmac251 frequencies measured between weeks 37-50, except for 133-06 which died 
early and data presented are from week 18. SIVsmE660 frequencies are measured between weeks 19-22. B) Phenotyping of cells from SIVmac251- 
(top, n = 4) and SIVsmE660-infected monkeys (bottom, n = 8) based on cell surface expression of CCR7 and CD28. Positivity for CCR7 and CD28 is 
indicated by the + and — signs below the bar graph. Cells were categorized as central memory (CCR7 + CD28 + ), transitional memory (CCR7~CD28 + ), or 
effector memory (CCR7~CD28~). All tetramer-positive cells were CD95 + . Phenotyping was conducted between weeks 44-78 for SIVmac251 and 33- 
46 for SIVsmE660. C) Measurement of percent of ex vivo perforin + granzyme B + cells from SIVmac251- (Left, n = 5) or SIVsmE660-infected monkeys 
(right, n = 8). D) Per-cell expression of perforin and granzyme B, measured by geometric mean of fluorescence (GMF) of perforin and granzyme B 
staining on cells from SIVmac251- (left, n = 5) and SIVsmE660-infected monkeys (right, n = 8). Perforin and granzyme B measurements were conducted 
between weeks 63-83 for SIVmac251 and 41-49 for SIVsmE660. In parts B-D, error bars represent mean ± SEM. * and **, significant at p<0.05 and 
0.01, respectively, after Bonferroni correction. NS, not significant. 
doi:10.1371/journal.ppat.1004069.g003 



technique developed in our laboratory that assesses the interaction 
of TCRs from a polyclonal epitope-specific CD8 T cell population 
with monomeric pMHC complexes [50] to evaluate the role of 
TCR:pMHC binding in determining the Mamu-A*0 1 -restricted 
epitope immunodominance hierarchy. 

For initial measurements of the binding of the pi 1C-, p54E660- 
, and p68A-specific TCRs to their respective monomeric pMHC 
complexes, we sorted total CD8 + T cells from seven Mamu-A *01 + 
SIVsmE660-infected monkeys. Detergent-resistant micro-domain 
(DRM) preparations from sorted cells, which were enriched in 



TCR complexes, were immobilized onto a Biacore LI chip. Initial 
assays demonstrated that specific binding to DRM preparations 
could be detected using monomers constructed with the dominant 
pi 1C and subdominant p54E660 epitope peptides (Figure 6A and 
Figure S8). The binding signal from the p54E660 monomers was 
weaker than that of the p 1 1 C monomers at the same concentra- 
tion. We were unable to detect specific binding of the monomers 
constructed with the more subdominant p68A epitope peptide at 
any of the concentrations of monomers (25-150 ug/mL) that were 
initially evaluated (Figure 6A). 
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Figure 4. The dominant pi 1 C-specific cells exhibited decreased antigen-specific expansion compared to subdominant epitope- 
specific cells. PBMCs from monkeys chronically-infected with either SIVmac251 (A) or SIVsmE660 (B) were stimulated in vitro with either p1 1C (red), 
p54E660/AS (blue), or p68A (green) peptide, harvested on days 3, 4, 5, 6, 8, 10, 12, and 14 following stimulation, and measured by flow cytometry to 
calculate the percent of tetramer-positive CD8 + T cells. Expansion was calculated as the fold change of the percent of each tetramer-positive 
population on each day, relative to day 0. Data from three SIVmac251- and three SIVsmE660-infected monkeys are shown. Measurements were 
conducted between weeks 40-52 for SIVmac251 and 31-44 for SIVsmE660. 
doi:10.1371/journal.ppat.1004069.g004 
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Figure 5. Dominant pi IC-specific population contained lower frequency of cytokine- and chemokine-producing cells than 
subdominant epitope-specific populations. PBMCs from monkeys chronically-infected with either SIVmac251 (n = 4) or SIVsmE660 (n = 5) were 
stimulated with either p1 1C, p54AS/E660, or p68A peptides and intracellular staining was used to assess production of the chemokine MIP-1 p and 
cytokines IFNy, TNFoc, and IL-2. Bars represent mean ± SEM. A) Individual cytokine and chemokine production. Top, SIVmac251. Bottom, SIVsmE660. 
B) Polyfunctional analysis. Positivity for each cytokine/chemokine is indicated by the dots below the bar graph. Vertical bars are grouped into 4, 3, 2 
or 1 function (indicated by the pink, light blue, purple, and orange horizontal bars, respectively). C) Polyfunctionality charts. Left, SIVmac251. Right, 
SIVsmE660. Each slice represents the percentage of tetramer-positive cells expressing between 1 and 4 functions. Data were collected between weeks 
36-42 for SIVmac251 and 14-25 for SIVsmE660. 
doi:1 0.1 371 /journal.ppat.1 004069.g005 



We performed titrations of the p 1 1 C and p54E660 pMHC 
monomers on DRM preparations from all seven monkeys 
(Figure 6B and Figure S9) and used a 1:1 Langmuir curve fitting 
analysis to calculate the association rate and the dissociation 
rate (kj) of the TCR:pMHC binding interaction (Table 3). The 
pi lC:Mamu-A*01 monomers bound to the DRM samples with a 
median k m of 7.38 x lOVMs (range 3.56-28.6 xl0 3 /Ms), which 
was faster than the p54E660:Mamu-A*01 monomers that had a 
median k m of 0.96xl0 3 /Ms (range 0.66-1. 87 xl0 3 /Ms). The 
dissociation rate of pi lC:Mamu-A*01 monomers from DRMs 
(median k g j- 0.024/s, range 0. 015-0. 030/s) was not substantially 
different from that of p54E660:Mamu-A*01 monomers (median 
kojf of 0.033/s, range 0.02 1—0.06 1 / s). The apparent dissociation 
constants (K d ) were derived from the kinetic rate constant values 
and showed that the pi lC:Mamu-A*01 monomers had a lower 
median K d (2.0 uM, range 1.0-6.8 uM) than the 
p54E660:Mamu-A*01 monomers (32 uM, range 22-43 uM). 
Thus, monomers constructed with the dominant p 1 1 C peptide 
epitope exhibited faster association rates to chip-immobilized 
DRMs, resulting in higher affinities (lower Kj) to TCRs than 
monomers constructed with the subdominant p54E660 epitope 
peptide. 

We speculated that the inability to detect the binding of 
p68A:Mamu-A*01 monomers to the DRM preparations may be 
due to a substantially lower affinity of the p68A:Mamu-A*0 1 
complex for its cognate TCRs. The difficulty in detecting binding 

Table 2. Expression of exhaustion-associated genes. 



Median Fold Change (p11C/p54AS) 



Gene (Protein) 


Day 14 


Day 70 


Week 31/32 


cd244 (2B4) 


-1.1 


-1.1 


1.0 


havcr2 (TIM-3) 


-1.0 


-1.1 


1.1 


Iag3 (LAG-3) 


3.2* 


1.0 


1.7* 


pdcdl (PD-1) 


-1.1 


-1.0 


-1.1 


ctla4 (CTLA-4) 


1.0 


1.0 


-1.0 


ctla4 


1.5* 


1.0 


-1.0 


ctla4 


-1.1 


-1.1 


1.0 


prdml (BLIMP-1) 


1.1 


1.1 


1.1 


prdml 


-1.0 


-1.2 


-1.0 


prdml 


1.0 


-1.1 


-1.0 


prdml 


-1.2 


-1.1 


-1.0 


cd160 (CD160) 


-1.1 


-1.2 


-1.2 


tbx21 (T-bet) 


1.1 


-1.1 


-1.1 


eomes (Eomesodermin) 


-1.0 


-1.7 


-1.1 


klrgl (KLRG1) 


-1.1 


1.1 


1.2 


klrg! 


1.1 


-1.1 


-1.0 



* pS0.05. 

doi:1 0.1 371 /journal.ppat.1 004069.t002 



was also likely exacerbated by the low frequency of the p68A- 
specific cells, which would have resulted in underrepresentation of 
the p68A-specific TCRs in the DRMs. Therefore, we pooled 
DRMs from p68A-specific cells sorted from multiple time points. 
These samples contained two- to six-fold more p68A-specific cells 
than samples from the initial experiments. 

Using these DRM preparations containing more p68A-specific 
cells and using higher concentrations of p68A:Mamu-A*0 1 
monomers (> 100 u.g/mL), we detected specific DRM binding 
(Figure 6C and Figure S10). However, specific monomer binding 
exhibited extremely fast association and dissociation rates and was 
detected only at very high concentrations. As a result, reliable K d 
values could not be measured using global curve fitting analyses. 
Our previous measurement of pllC and p54E660 pMHC 
monomer binding to DRMs indicated that p 1 1 C had the fastest 
association rate. The association rate of the p68A monomer to its 
respective TCRs was at least as fast as the association of the 
quantifiable p 1 1 C monomer, although likely to be even faster. 
Therefore, we estimated that the k on for p68A was at least 10x10 / 
Ms. Knowing that the Biacore instrument's lower limit of 
detection of dissociation is about 1 /s, we estimated that the k a ff 
of the p68A:Mamu-A*0 1 monomer from its respective TCRs was 
at least as fast as 1 /s. Using these estimated k on and k 0 j/-, the 
apparent of the p68A monomer was estimated be at least 
100 uM (Table 3). 

Therefore, the calculated K d values were lowest for pllC 
monomers, intermediate for p54E660 monomers, and estimated 
to be highest for p68A monomers. In summary, epitope 
dominance was associated with the hierarchy of TCR affinities. 
Affinity differences were driven by the fast association rate of the 
dominant pllC epitope:Mamu-A*01 complex to its respective 
TCRs and likely fast dissociation rate of the more subdominant 
p68A:Mamu-A*0 1 complex. 

Discussion 

The development of an effective vaccine against a pathogen 
such as HIV-1 will require a detailed understanding of the 
interplay between epitope immunodominance and the functional 
capacities of epitope-specific CD8 + T cells. The present study was 
undertaken to explore this relationship in SIV-infected rhesus 
monkeys. 

It is important to note that our longitudinal gene expression 
analyses were conducted with a restricted sample size of six 
monkeys. In addition, because of the large numbers of cells 
required for RNA isolation, we were unable to perform 
phenotyping or functional assays during acute infection. Despite 
these clear limitations, we were able to detect a number of 
significant differences in gene expression between the dominant 
and subdominant epitope-specific populations during acute 
infection. Baron et al had similar findings, in a murine model, 
where they characterized a dominant (H7' 1 ) and a cryptic (HY) 
H2-D b -restricted epitope-specific CD8 + T cell population [20]. In 
the murine model, the cryptic epitope-specific cells had higher 
IL7R expression during primary expansion and higher L-selectin/ 
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Figure 6. TCR:pMHC measurements. DRMs were purified from total CD8 + T cells sorted from seven chronically-infected SIVsmE660-infected 
monkeys. The DRMs were evaluated for specific binding, measured in resonance units (RU), to pMHC monomers constructed with pi 1 C, p54E660, 
and p68A epitope peptides and Mamu-A*01. Data are representative of the binding observed from all seven monkeys evaluated. A) Initial 
experiments to detect p1 1C, p54E660, and p68A monomer binding to DRMs. 100 ug/mL p1 1C (red), 100 u.g/ml_ p54E660 (blue), and 150 ug/mL p68A 
(green) pMHC monomer binding are overlaid from experiments conducted on separate Biacore Chips. Binding of the p68A:Mamu-A*01 monomer to 
the DRMs was not observed at any concentration of monomer tested (25-200 ug/mL). B) Titrations of pi 1 C (top) and p54E660 (bottom) 
peptide:Mamu-A*01 monomers for calculation of binding kinetics and affinity. Overlaid sensograms of the binding of pi 1 C and p54E660 pMHC 
monomers to DRMs purified from total CD8 + T cells are shown. Monomer binding was evaluated using pMHC concentrations ranging from 25 to 
200 ug/mL The black curve shows the Langmuir fitted curve that was used to calculate kinetics. C) Detection of p68A peptide:Mamu-A*01 monomer 
binding to DRMs. p68A-specific CD8 + T cells were collected from multiple tetramer-specific flow cytometric cell sorts and pooled for DRM purification. 
Titrations of p68A pMHC monomers were performed at concentrations ranging from 150 to 1000 ug/mL Binding to DRMs from monkey ZD57 at 
1000 ug/mL is shown and is representative of the four monkeys that were evaluated. All readings have been normalized by subtracting the binding 
of the control monomer TL8 run at the same concentrations as the experimental monomers. 
doi:1 0.1 371 /journal.ppat.1 004069.g006 



SELL on all timepoints evaluated. Similarly, we found that both 
IL7R and SELL were higher in the subdominant p54AS-specific 
cells. Additionally, Baron et al found that G^MA was highly 
expressed in the dominant epitope-specific CD8 + T cells. We 
found a trend toward higher expression of GJ(MA and significandy 
higher expression of G^AIB and G^MK in the dominant epitope- 
specific populations. This suggests that common gene networks are 
involved in the regulation of dominant immunologic responses. 

We also evaluated the expression of exhaustion-associated genes 
between the dominant and subdominant epitope-specific cells. We 



found that the dominant pi lC-specific cells had increased 
expression of genes encoding LAG-3 and CTLA-4 during acute 
infection, both of which are known to be upregulated upon TCR 
stimulation [51,52]. The role of these exhaustion-associated genes 
during acute infection is not well established and it is unclear to 
what extent these molecules impair CD8 + T cell function during 
acute HIV/SIV infection. Their increased expression on the 
dominant pi lC-specific cells during acute infection suggests that 
these cells were experiencing a greater level of activation driven by 
greater stimulation through their higher-affinity TCRs. On day 
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Table 3. TCR:pMHC binding values. 







*„„ (xio'wrV 1 




Kd (filVl) 


plIC 


7.38 (3.56-28.6) 


0.024 (0.015-0.030) 


2.0 (1.0-6.8) 


p54E660 


0.96 (0.66-1.87) 


0.033 (0.021-0.061) 


32 (22-43) 


p68A 


aio* 


atl* 


2:100* 



Values are the median (range) from seven monkeys. 
k on , association rate. 
k off , dissociation rate. 

K d , equilibrium dissociation constant. Calculated using the equation k off /k on . 
* k on of the p68A monomers was estimated as being at least as fast as the 
fastest measured monomer, p11C. 

*k off of p68A monomers was estimated as being at least as fast as the limit of 
detection of the Biacore instrument, approx.l/s. 

1 K d of p68A monomer binding calculated using estimated values of the k on and 

k off using the equation k off lk 0 „. 

doi:10.1371/journal.ppat.1004069.t003 

70, no differences in gene expression were noted between the 
dominant and subdominant epitope-specific cells. However, by 
late chronic infection, the dominant p 1 1 C-specific cells showed 
increased levels of LAGS transcripts. LAG-3 is a well-described 
marker of exhausted CD8 + T cells during chronic viral infections 
[53,54], including HIV [55]. However, while LAG-3 expression 
on CD8 + T cells has been shown to contribute to immune 
suppression in certain mouse models [56], its suppressive effect 
during HIV/SIV infection remains unclear [57,58]. Additional 
studies are needed to confirm the cell surface expression of this 
protein and determine the extent to which increased expression 
contributes to functional differences between dominant and 
subdominant S1V/H1V epitope-specific cells. 

In the present study, we found that dominant epitope-specific 
cells had greater cytotoxic potential, whereas the subdominant 
epitope-specific cells displayed enhanced proliferation capacity 
combined with increased cytokine and chemokine production. 
Importantly, these dominant epitope-specific cells had higher 
TCR affinity, a measure that typically correlates with functional 
avidity [59-61]. The ability of HIV/SIV epitopes with higher 
functional avidity to induce T cells that are less polyfunctional with 
lower proliferative capacity has been described by Harari et al [62] 
and Vigand et al [63] . These studies also corroborate our findings 
that cells of increased avidity also exhibited increased expression of 
exhaustion-associated molecules. Similarly, Conrad et al found that 
dominant TCR clonotypes exhibited reduced in vitro cytokine 
production and survival compared to the subdominant clonotypes 
that was associated with increased TCR avidity and increased 
expression of PD-1 [64]. These studies support the notion that 
reduced CD8 + T cell avidity promotes the development of 
polyfunctional T cells that under conditions of persistent high 
antigenemia may be able to partially control viral replication. This 
is also supported by the fact that increased avidity is significandy 
associated with viral escape [65-67] and increased induction of 
cell death and exhaustion [68,69]. Collectively, these studies are 
consistent with our findings and suggest that to generate 
polyfunctional T cells after vaccination, the antigens used cannot 
confer extremely high TCR affinity; otherwise exhausted and 
monofunctional T cells will be generated. 

The finding that the dominant pllC-specific cells were not 
highly polyfunctional is surprising considering that these cells are 
thought to be predominantly responsible for control of SIV viral 
replication [70,71]. This reduced polyfunctionality may be 
compensated by a number of other cell characteristics. First, they 
may possess higher functional avidity, a characteristic that also 



correlates strongly with HIV control [15,60,72,73] and higher 
TCR affinity [59-61]. Second, the dominant epitope-specific cells 
are present at substantially higher frequencies and contain greater 
amounts of perform and granzyme than the subdominant epitope- 
specific cells. Thus on a percentage basis, the former population 
has an increased capacity to kill virus-infected cells. Finally, they 
may have greater TCR cross-reactivity to epitope variants, a 
characteristic that is thought to be an important component of 
durable HIV control by CD8 + T cells [73,74]. In support of this 
possibility, Mothe et al found that HIV controllers had higher 
functional avidity that was associated with higher cross-reactivity 
of CD8 + T cells to HIV epitope variants compared to non- 
controllers [73]. Additionally, Bennett et al generated a panel of 
mutant TCRs specific for the HIV epitope SL9 and found that 
those TCR mutants with higher TCR avidities generally were 
more cross-reactive to SL9 variants than those with lower avidities 
[61]. In contrast, the previously-mentioned study by Conrad et al 
found that the high TCR avidity of the dominant clonotypes had 
reduced cross-reactivity than subdominant clonotypes [64] . Future 
studies should compare the TCR cross-reactivity to variant 
epitopes between the dominant plIC- and the subdominant 
epitope-specific cells to determine the extent that it contributes to 
the protective capacity of the pi 1 C-specific cells. 

We found that TCR affinity was predictive of frequency of the 
dominant and subdominant epitope-specific cells that are gener- 
ated. The role of TCR:pMHC interactions in determining 
immunodominance hierarchies has been largely unexplored. A 
limited number of murine studies that have addressed this by 
measuring TCR dissociation rates using peptide:MHC class I 
tetramers have either found a lack of association between TCR 
affinity and epitope dominance [75] or found an association in the 
opposite direction [76]. The discordance between these studies 
and ours may be that for the measurement of the TCR 
dissociation rate the mouse studies used pMHC tetramers as 
opposed to monomers used in the present study. It is known that 
the accuracy of pMHC tetramers to predict physiologic 
TCR:pMHC binding affinities is limited [77,78]; The multi- 
valency of tetramers combined with their propensity to aggregate 
can complicate accurate measurement of these interactions [79]. 
Moreover, the present study found that overall binding affinity 
(which incorporates both association and dissociation rate), rather 
than dissociation rate alone, was predictive of epitope dominance. 
In the above-mentioned murine studies, only the dissociation rate 
was assessed, leaving the association rate and affinity unknown. 
Slower association rates have been previously attributed to 
induced-fit or thermodynamic entropic penalty [80-82] in TCR 
binding and as such, pMHC epitopes that bind with faster 
association rates may outcompete those that bind with slower 
kinetic rates and gain dominance in response. Therefore, it is 
possible that association rate also played a role in the 
immunodominance hierarchies in these mouse studies and its 
measurement may have been informative about the particular 
parameter of the TCR:pMHC interactions that was important for 
immunodominance. 

The present study found that high-frequency dominant epitope- 
specific CD8 + T cells had reduced proliferation and cytokine 
production but increased cytotoxicity potential compared to low- 
frequency subdominant epitope-specific CD8 + T cells. This was 
associated with higher TCR affinity and increased expression of 
exhaustion-associated genes in the dominant epitope-specific cells. 
The engineering of higher-affinity TCRs has been explored for 
enhancement of HIV immunity [61,83]. However, our findings 
suggest that high TCR affinity may increase susceptibility to 
exhaustion and reduced functionality, although it may also 
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promote higher cytotoxic potential. Since both proliferation and 
cytokine production as well as cytotoxicity are associated with 
superior control of HIV replication, antigen targets of engineered 
TCRs or antigens that are included in an HIV vaccine should be 
of both high and low affinity. The high-affinity antigens will 
promote epitope-specific CD8 + T cells of high frequency with high 
cytotoxic potential while the low-affinity antigens will promote 
CD8 + T cells that are more polyfunctional and proliferative and 
that have decreased potential for exhaustion. 

Materials and Methods 

Ethics statement 

Indian-origin rhesus monkeys used in this study were main- 
tained according to the guidelines of the National Institutes of 
Health (NIH) Guide to the Care and Use of Laboratory Animals 
and the approval of the Institutional Animal Care and Use 
Committee (IACUC) of Harvard Medical School (protocol # 
03503) and the NIH. The institution also accepts as mandatory the 
Public Health Service Policy on Humane care and use of 
Laboratory Animals by Awardee Institutions and the NIH 
Principles for the Utilization and Care of Vertebrate Animals 
used in Testing, Research and Training. The New England 
Primate Research Center (NEPRC) is fully accredited by the 
Association for Assessment and Accreditation of Laboratory 
Animal Care International (AAALAC). The NEPRC has devel- 
oped a comprehensive environmental enrichment, psychological 
well-being, and minimization of distress plans for primates that is 
available for inspection by the United States Animal and Plant 
Health Inspection Service (APHIS) and to officials of any pertinent 
funding agency. The Harvard Medical School IACUC documents 
NEPRC compliance with the plan during semiannual facility 
inspections. 

Animals 

Mamu-A*01-po$itive, Mamu-B*l 7-nega.tive, Mamu-B*08-nega- 
tive, Indian-origin rhesus monkeys (Macaca mulatto) were selected 
for these studies by PCR-based MHC typing as previously 
described [84] . All monkeys were infected intrarectally with either 
SIVmac.251 or SIVsmE660. SIVmac251 infection was adminis- 
tered by a single high dose challenge. SIVsmE660 was adminis- 
tered by a repeated low-dose challenge. 

Antibodies and flow cytometry reagents 

Conjugated antibodies and staining reagents included MlP-lfi- 
PE, CD3-PB, CD3-PE-Cy7, CD3-PerCP-Cy5.5, CD3-APC-Cy7, 
CD3-Horizon V450, CD4-PerCP-Cy5.5, CD4-AmCyan, CD4- 
FITC, CD8a-APC, CD8a-APC-Cy7, CD8a-AlexaFluor700, 
CD8a-FLTC, CD8ot-APC-H7, CD69- ECD (Beckman Coulter), 
CD20-Horizon V450, CCR7-FITC (R&D Systems), CD95-APC, 
CCR7-PerCP-Cy5.5, CD28-PE-Cy7 (eBioscience), granzyme B- 
AlexaFluor700, perforin-FITC (MabTech), IFNy-PE-Cy7, TNFot- 
AlexaFluor700, IL-2-APC, CD95-PE, CD95-APC, and Aqua 
LIVE/DEAD Fixable Dead Cell Stain (Invitrogen). All reagents are 
from BD Biosciences unless indicated otherwise. For construction of 
monomers and tetramers, the following peptides were synthesized 
and purified to >95% by HPLC by New England Peptide LLC: 
pllC (CTPYDINQM), p54AS (TVPWPNASL), p54E660 
(TVPWPNETL), p68A (STPPLVRLV), and TL8 (TTPESANL). 
The monomers and tetramers were prepared as previously 
described [85,86]. Tetramers were prepared using either strepta- 
vidin-PE (Prozyme), -APC (Prozyme), -AlexaFluor488 (Invitrogen), 
or -Qdot655 (Invitrogen). Monomers used in surface plasmon 



resonance studies were further quantified using an RC DC protein 
kit (Bio-Rad). 

RNA extraction and microarray analysis 

Peripheral blood was collected on days 14, 21, 28, 35, 42, 56, 
70, and at week 31 or 32 post-inoculation with SIVmac251. 
Plasma viral RNA from these samples was measured using an 
ultra-sensitive branched DNA amplification assay (Siemens 
Diagnostics, Berkeley, CA). PBMC were stained with pllC and 
p54AS tetramers and CD3 and CDS antibodies at 4°C. Tetramer- 
positive single CD3 + CD8 + lymphocytes were sorted to &95% 
purity into RNAprotect (Qiagen) at 4°C. RNA was isolated using a 
Trizol (Invitrogen) extraction protocol. Briefly, 0.8 ml of Trizol 
were added to the cell pellet and incubated for 5 minutes at room 
temperature and 0.16 ml of chloroform were added, shaken 
vigorously by hand for 1 5 seconds, incubated at room temperature 
for 2-3 minutes and centrifuged at 13,000 rpm for 15 minutes at 
4°C. The colorless upper aqueous phase was carefully collected 
and transferred to a new tube containing 2 ul of linear acrylamide 
for mixing. An equal volume of isopropyl alcohol was then added 
and mixed. The mixture was incubated at room temperature for 
10 minutes and centrifuged at 13,000 rpm for 10 minutes at 4°C. 
The supernatant was collected and the RNA was washed with 
1 ml of 70% ethanol and centrifuged at 10,500 rpm for 5 minutes 
at 4°C. The supernatant was completely removed and the RNA 
pellet was allowed to air-dry. The RNA was then resuspended in 
RNase-free water and stored at — 80°C. RNA integrity was tested 
using an Agilent Bioanalyzer. RNA was then amplified using the 
TargetAmp 2-Round Biotin-aRNA Amplification Kit 3.0 (Epi- 
centre Biotechnologies) according the manufacturer's instructions. 
Amplified biotinylated antisense-RNA (aRNA) was resuspended in 
RNase-free water and stored at — 80°C. Nanodrop ND-1000 was 
used to determine the biotinylated aRNA concentration and an 
Agilent Bioanalyzer was used to determine its integrity. Amplified 
aRNA was hybridized to Illumina Human HT-12 Expression 
BeadChips according to the manufacturer's instructions and was 
stained with Streptavidin Cy3 for detection (Illumina, San Diego, 
CA, USA). The Human HT-12 BeadChip assays 48,000 
transcripts. The BeadChips were built with sequences derived 
from the National Center for Biotechnology Information Refer- 
ence Sequence (NCBI RefSeq) database (Build 36.2, Release 22). 
Arrays were scanned according to the manufacturer's instructions. 
Processing of the raw array data was performed using Illumina 
BeadStudio software. Matlab (Mathworks, Natick, MA, USA) was 
used to perform statistical analysis of gene expression data. The 
entire expression dataset was first quantile-normalized. Fold- 
differences of individual transcripts were computed as the ratios of 
the median expression levels of the pi 1C- to those of the p54AS- 
specific CD8 + T cells (pi 1C/ P 54AS). When p54AS expression 
values were larger than pllC values, the negative reciprocal was 
calculated. A Wilcoxon signed-rank test was used to evaluate the 
significance of the differences in gene expression between pllC- 
and p54AS-specific CD8 + T cells. Genes whose expression was 
significandy different between the pllC- and p54AS-specific cell 
samples and whose median fold change difference was at least 
±1.5 were considered differentially expressed. The Minimum 
Information about a Microarray Experiment (MIAME) criteria 
have been met by these experiments. The complete dataset is 
available in the Gene Expression Omnibus, accession number 
GSE54001. 

Flow cytometric analysis 

In order to collect and analyze comparable numbers of each 
CD8 + T cell specificity, different numbers of PBMCs were used for 
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each pllC, p54, and p68A sample. Input numbers of PBMCs 
were dependent on the expected relative frequency of each 
epitope-specific population. For example, if the expected frequen- 
cies of the pllC-, p54-, and p68A-specific CD8 T cells were 
10%, 1%, and 0.1%, respectively, then the ratio of the numbers of 
cells used for each sample was 1:10:100. Staining volumes, 
washing volumes, and amounts of staining reagents used were 
adjusted for each sample to ensure the same cell-to-reagent 
concentration ratios were used for all samples. For quantification 
of epitope-specific CD8 + T cells, PBMCs were stained with 
tetramer and CD3, CD4, and CDS antibodies. For phenotyping, 
PBMCs were additionaUy stained with CD20, CD95, CD28, and 
CCR7 antibodies. For ex vivo measurement of perforin and 
granzyme B, PBMCs were stained with tetramer, CD3, CD4, and 
CDS antibodies and then stained intracellularly with perforin and 
granzyme B antibodies as previously described [87]. FMO samples 
were also prepared to set the gates for perforin and granzyme B 
gating. All flow cytometric data were collected on an LSRII (BD) 
and analyzed using Flowjo (Tree star). Statistical analysis of flow 
cytometric data was conducted in GraphPad Prism 5. Wilcoxon 
signed rank test was used to compare groups. P-values SO. 05 were 
considered significant. Bonferroni correction was used for multiple 
comparisons with significant p-values ^0.05/(number of compar- 
isons) considered significant. 

Cytokine and chemokine production 

Measurements of IL-2, IFNy, TNFoc, and MIP-ip production 
by epitope-specific CD8 T cells following stimulation of 1 nM 
peptide was performed in 5 mL FACS tubes as previously 
described [50]. Flow cytometric data were gated on single, live, 
CD3 + CD4~CD8 + CD69 + lymphocytes. Boolean gating was used 
to calculate percent of polyfunctional cells. 

Expansion of epitope-specific CD8 + T cells 

PBMCs, at approximately 3xl0 6 lymphocytes/mL (measured 
using Guava EasyCyte automatic cytometer (Millipore)), were 
plated in flat-bottomed culture plates and stimulated for expansion 
with 1 nM epitope peptide. Approximately lxlO 6 PBMCs were 
also stained with tetramers and CD3, CD4, and CDS antibodies, to 
measure the frequency each CD8 + T cell specificity before peptide- 
stimulated expansion (day 0). 20 U/mL IL-2 (Hoffmann- LaRoche) 
was added to stimulated cells on day 3, and media supplemented 
with IL-2 was changed thereafter as needed. Samples were 
harvested on days 3, 4, 5, 6, 8, 10, 12, and 14 post-stimulation for 
staining with aqua LIVE/DEAD, tetramer and CD3, CD4, and 
CDS antibodies for flow cytometric analysis. Flow cytometric data 
were gated on single, live, CD3 + CD4 CD8 + tet + lymphocytes. 

Surface plasmon resonance (SPR) measurements 

Freshly isolated PBMCs were enriched for CD8 + T lymphocytes 
using Miltenyi's magnetic-activated cell sorting (MACS) kit for 
isolation of untouched NHP CD8 + T cells according to the 
manufacturer's instructions. Isolated CD8 + T cells were then 
stained at 4°C with tetramers and CD3 and CDS antibodies. 
Sorting was performed on an Aria flow cytometer/cell sorter (BD). 
For initial TCR binding analyses, cells were only stained with the 
TL8 tetramers and those cells that were positive for TL8 were 
excluded from the sorted CD3 + CD8 + lymphocytes. For repeated 
TCR binding analysis for enrichment of p68A-specific TCRs, cells 
were stained with TL8, p54E660, and p68A tetramers, and 
multiple populations were simultaneously sorted: p54E660 
p68A" TL8", p54E660 + p68A" TL8", and p54E660" p68A + 
TL8 . These sorted populations were subsequently combined 
during purification of the DRMs. For samples with low total cell 



number, extra DRMs from cells not specific for any of the epitopes 
being evaluated were added to provide extra mass to the pellet 
during centrifugation. These extra DRMs were obtained from 
sorted single CD3 CDS lymphocytes that were negative for all 
evaluated epitope specificities (pi lC~p54E660~p68A~TL8~). All 
sorted cells were resuspended in a 4°C solution of protease 
inhibitors (1 pg/mL each leupeptin, pepstatin, and aprotinin) and 
stored at — 80°C until analysis. DRM purification was carried out 
as described previously [50] with the modification that the repeat 
samples for measurements of p68A were conducted at concentra- 
tions up to 1,000 pg/mL of monomer. The specific binding signal 
was obtained by subtracting the non-specific signal from TL8 
pMHC control monomer binding to the TL8-depleted DRM 
preparation when injected at the same concentration from the 
pllC, p54E660, or p68A binding signal. The global curve fitting 
to the Langmuir equation was used to derive kinetic rate constants 
(k m and k 0 Jj for calculation of equilibrium dissociation constant Kj 
as previously described [88,89]. All SPR measurements were 
carried out on a Biacore 3000 instrument, and data analyses were 
performed using BIAevaluation 4.1 software (GE Healthcare). 

Supporting Information 

Figure SI Gating of pllC- and p54AS-specific CD8 + T 
cells sorted for microarray analyses. Gates for pllC (y- 
axis) and p54AS (x-axis) tetramers on samples sorted from each of 
the six Mamu-A*0t rhesus monkeys (107-06, 112-06, 125-06, 133- 
06, 142-06, and 202-06) used in the gene expression analysis. Plots 
are gated on single CD3 + CD8 + lymphocytes and percent of 
tetramer-positive cells are shown on each plot for the p 1 1 C- (upper 
left) and p54AS- (middle right) positive cells. 50,000 CD8 + T cells 
events are displayed for all samples except where indicated in 
upper right corner of plots. The data file for 107-06 on day 56 was 
corrupt and not available for re-analysis. 
(TIF) 

Figure S2 Individual values of genes differentially 
expressed between dominant pllC- and subdominant 
p54AS-specific CD8 T cells. Shown are the individual 
normalized expression, measured in fluorescence units, for each 
A) maturation, B) cytotoxicity, and C) proliferation and apoptosis 
gene that was determined to be differentially expressed between 
pllC- and p54AS-specific cells. Plots include values for which its 
matching pair is missing and therefore were not used in 
determination of differential expression. Fold change values are 
indicated in upper left corner. An asterisk under fold change values 
is present if that gene met the criteria for differential expression on 
that timepoint. For those genes with more than one probe on the 
BeadChip, the numerical probe IDs are included in the gene 
name. Red, pi 1C. Blue, p54AS. 
(PDF) 

Figure S3 Staining and gating of tetramers for quanti- 
fication of epitope-specific cells in Mamu-A*01 immu- 
nodominance hierarchies. Gating of pllC-, p54AS/E660-, 
and p68A-specific CD8 + T cells from A) SIVmac251- and B) 
SIVsmE660-infected animals. 50,000 CD8 + T cell events are 
displayed on each flow plot. SIVmac251 frequencies were 
measured between weeks 37-50, except for 133-06 which died 
early and data used are from week 18. SIVsmE660 frequencies 
were measured between weeks 19-22. 
(TIF) 

Figure S4 Representative flow staining and gating for 
phenotyping. Representative staining and gating for phenotyp- 
ing of tetramer-positive cells from (A) SIVmac251- and (B) 
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SIVsmE660-infected monkeys. Left, Staining and gating of 
tetramer-positive cells among CD8 + T cells (40,000 CD8 + T cell 
events are shown). Right, staining and gating of CD28 (y-axis) and 
CCR7 (x-axis) cell surface proteins among the tetramer-positive 
cells (1,000 tetramer-positive events are shown). Phenotyping was 
conducted between weeks 44-78 for SIVmac251 and 33-46 for 
SIVsmE660. 
(TIF) 

Figure S5 Representative staining and gating of per- 
forin and granzyme B. Mamu-A*01 tetramers were used to 
identify pllC-, p54AS/E660-, and p68A-specific cells among 
PBMCs from monkeys chronically-infected with either SIV- 
mac251 or SIVsmE660. Perforin and granzyme B content within 
the tetramer-positive cells were measured by flow cytometry. A) 
FMOs used to set gates for perforin (left) and granzyme B (right). 
For SIVmac251 (B) and SIVsmE660 (D), representative staining 
and gating of tetramer-positive cells (left) and perforin + granzyme 
B + cells among the tetramer-positive cells (right). 250,000 CD8 + T 
cells are shown on flow plots of tetramer staining. 1,000 tetramer- 
positive events are shown on flow plots of perforin and granzyme B 
staining. For SIVmac251 (C) and SIVsmE660 (E), histograms 
illustrating fluorescence intensity of perforin (top) and granzyme B 
(bottom) within the pllC (red), p54AS/E660 (blue), and p68A 
(green) tetramer-positive cells. FMOs shown as filled gray 
histograms. Measurements were conducted between weeks 63- 
83 for SIVmac251 and 41-49 for SIVsmE660. 
(TIF) 

Figure S6 Representative staining and gating of tetra- 
mers during in vitro expansion. PBMCs from monkeys 
chronically-infected with either SIVmac25 1 (A) or SIVsmE660 (B) 
were stimulated in vitro with either pllC, 54E660/AS, or p68A 
peptide, harvested on days 3, 4, 5, 6, 8, 10, 12, and 14 following 
stimulation, and measured by flow cytometry to calculate the 
percent (upper right on each flow plot) of tetramer-positive CD8 + 
T cells. Expansion was calculated as the fold change (middle left 
on each flow plot) of the percent of each tetramer-positive 
population on each day, relative to day 0. Number of CD8 + T cell 
events displayed are the same for pllC, p54AS/E660, and p68A 
plots within a given day. Measurements were conducted between 
weeks 40-52 for SIVmac251 and 31-44 for SIVsmE660. 
(TIF) 

Figure S7 Representative staining and gating for cyto- 
kine and chemokine production. PBMCs from monkeys 
chronically-infected with either SIVmac25 1 (A) or SIVsmE660 (B) 
were stimulated with either pi 1C, p54AS/E660, or p68A peptides 
and intracellular staining was used to assess production of the 
chemokine MIP-10 and the cytokines IFNy, TNFot, and IL-2. 
Left, representative staining and gating of tetramer-positive cells 
(25,000 CD8 + T cell events are shown). Right, corresponding 
staining and gating of MIP-ip, IFNy, TNFot, and IL-2 (250 
tetramer-positive events are shown). Individual gates were then 
used in a Boolean analysis for assessment of polyfunctionality. 
Data was collected between weeks 36-42 for SIVmac251 and 14— 
25 for SIVsmE660. 
(PDF) 
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Figure S8 Detection of specific binding of pllC and 
p54E660 peptide:Mamu-A*01 monomers to DRMs. 

DRMs were purified from total CD8 + T cells sorted from seven 
chronically-infected SIVsmE660-infected monkeys. The DRMs 
were evaluated for specific binding, measured in resonance units 
(RU), to pMHC monomers constructed with pi 1C, p54E660, and 
p68A epitope peptides and Mamu-A*01. Shown are overlaid 
readings of the binding of pi 1C (red) and p54E660 (blue) pMHC 
monomers at 100 ug/mL. p68A:Mamu-A*0 1 monomer binding 
above background was not detected at any concentration and is 
not shown. Readings have been normalized by subtracting the 
binding of the control monomer TL8 run at the same 
concentrations as the experimental monomers. 
(TIF) 

Figure S9 Titrations of pllC and p54E660 peptide:- 
Mamu-A*0 1 monomers for calculation of binding kinet- 
ics and affinity. Shown are sensograms indicating the binding 
of pllC (left) and p54E660 (right) pMHC monomers to DRMs 
purified from total CD8 + T cells sorted from seven SIVsmE660- 
infected monkeys. pllC monomers were run at 25 (green), 50 
(pink), 100 (blue), and 200 (red) ug/mL. The ARB0 plot for pi 1C 
shows a 150 |ig/mL (light blue) run in place of the 100 ug/mL. 
The AP34, ZD57, and A6V031 plots for pi 1C do not show the 
200 ug/mL run. p54E660 monomers were run at 25, 50, 150, and 
200 ug/mL for AP54, ARB0, 8B1, and AS47 and at 25, 100, and 
200 ug/mL for AP34, ZD57, and A6V031. The ZD57 plot 
includes an additional 50 |tg/ mL run. A Langmuir curve was fit to 
each binding curve at each concentration and was used to 
calculate binding kinetics. Readings have been normalized by 
subtracting the binding of the control monomer TL8 run at the 
same concentrations as the experimental monomers. 
(PDF) 

Figure S10 Detection of p68A monomer binding. p68A- 
specific CD8 T cells were collected from multiple tetramer- 
specific flow cytometric cell sorts and pooled for DRM 
purification. Titrations of p68A pMHC monomers were per- 
formed at concentrations ranging from 150 to 1000 ug/mL. The 
highest concentration evaluated is shown. Binding of the control 
monomer TL8 at the same concentration has been subtracted 
from all readings. 
(TIF) 
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